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Abstract. We present results derived from VLT-FORS2 spectra of 24 different globular clusters associated with 
the lenticular galaxy NGC 3115. A subsample of 17 globular clusters have sufficiently high signal-to-noise to allow 
precision measurements of absorption line-strengths. Comparing these indices to new stellar population models 
by Thomas et al. we determine ages, metallicities and element abundance ratios. For the first time these stellar 
population models explicitly take abundance ratio biases in the Lick/IDS stellar library into account. Our data 
are also compared with the Lick/IDS observations of Milky Way and M31 globular clusters. Unpublished higher 
order Balmer lines (H7a,f and H5a,f) from the Lick/IDS observations are given in the Appendix. 
Our best age estimates show that the observed clusters which sample the bimodal colour distribution of NGC 3115 
are coeval within our observational errors (2-3 Gyr). Our best calibrated age/metallicity diagnostic diagram (H/3 
vs [MgFe]) indicates an absolute age of 11-12 Gyr consistent with the luminosity weighted age for the central 
part of NGC 3115. We confirm with our accurate line-strength measurements that the {V — I) colour is a good 
metallicity indicator within the probed metallicity range (—1.5 < [Fe/H] < 0.0). The abundance ratios for globular 
clusters in NGC 3115 give an inhomogeneous picture. We find a range from solar to super-solar ratios for both 
blue and red clusters. This is similar to the data for M 31 while the Milky Way seems to harbour clusters which are 
mainly consistent with [a/Fe] ~ 0.3. From our accurate recession velocities we detect, independent of metallicity, 
clear rotation in the sample of globular clusters. 

In order to explain the metallicity and abundance ratio pattern, particularly the range in abundance ratios for the 
metal rich globular clusters in NGC 3115, we favour a formation picture with more than two distinct formation 
episodes. 

Key words, galaxies: abundances - galaxies: individual (NGC 3115) - galaxies: star clusters - galaxies: stellar 
content - globular clusters: general 



1. Introduction 



lapse and merger models for galaxy formation (e.g., Forbes 



The analysis of globular cluster systems in external galax- 
ies is starting to fulfil its long-held promise as a probe 
of the formation of galaxies. Globular clusters (hereafter 
GCs) are a fossil record of this formation process, and 
provide one of the best tools with which to investigate 
the chemical enrichment and star formation history in the 
initial stages of galaxy formation (e.g.. Ashman & Zepf 
1998| ). Much recent interest has focused on the globular 
cluster systems of luminous elliptical galaxies where the 
combination of metallicities and kinematics can be used to 
distinguish between variants of the popular monolithic col- 



et al.||l997l ; [Sharpies et al.||1998| ) 



* Based on observations collected at the European Southern 
Observatory, Cerro Paranal, Chile (ESO No. 66.B-0131). 



Lenticular galaxies hold a key position in the Hubble 
sequence of morphological types, intermediate between 
pure spheroidal systems like luminous ellipticals and disk- 
dominated spiral galaxies. Their formation mechanism is 
still the subject of considerable debate with evidence both 
for (Dressier et al. 1997) and against (Dressier 198C| ) their 
evolution from star-forming spirals via processes of gas 
stripping and exhaustion. A key question is when and how 
did such processes occur for SO galaxies in a wide range 
of environments from rich clusters to the field. The globu- 
lar cluster systems of SO galaxies can provide independent 
constraints on when the major star formation episodes oc- 
curred both in the disk and halo. However, thus far they 



have been little studied with only NGC 1380 (Kissler-Patig 
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et al. |1997|) and NGC 4594 (the Sombrero galaxy, [Bridges 
et al. |L997 ) having received any detailed attention and 
only with photometric methods. 

NGC 3115 is one of the nearest SO galaxies (9.7 ± 
0.4 Mpc, Mb -20.1; iTonry et al.l |200l|) and is lo- 



cated in the sparse low-density environment of the Leo 
Spur. As such it provides an ideal test case for study- 
ing the formation mechanism of field SO's. A significant 
globular cluster system containing ^^500 clusters was first 
detected by Hanes & Harris (198(;) using photographic 



plates. The nature of the cluster system and its origin 
were recently thrown into question with the discovery by 



Elson] ( |1997D that the red-giant stars in the NGC 3115 halo 
~40 kpc from the centre showed a himodal colour distri- 
bution. The inferred presence of two distinct halo popula- 
tions of roughly equal size at metallicities of [Fe/H] ~ —0.7 
and [Fe/H] ~ — 1.3 suggests at least two distinct epochs 
of formation. 

The {V - 1) colour distribution of the NGC 3115 glob- 
ular cluster system has been the target of two recent in- 
dependent st udies using HST ( Kundu fc Whitmorc 1998 ) 
and CFHT ( [Kavelaar^ |199^ ) data. Both studies find bi- 
modality in the colour distributions of the GCs, with mean 
metallicities at [Fe/H] ~ —0.37 and [Fe/H] ~ —1.36 sug- 
gesting that the cluster and halo star systems may have 
formed coevally. This suggestion has gained further sup- 
port from an investigation by Puzia et al. (2002b) who 



employed optical-IR colours to probe the globular clus- 
ter population close to the centre of NGC 3115. They also 
find two peaks in metallicity and an average age around 
~10 Gyr. However, their age discrimination power is very 
limited for metallicities lower than [Fe/H] = —0.4. 

One scenario in which the above observations could be 
understood is if the metal-poor component corresponds to 
a primordial ~ 13 Gyr old population, whilst the metal- 
rich component formed a few Gyr later from enriched gas , 
possibly as the result of a minor merger (e.g., Bekki 1998| ). 
With only broad-band colours available, however, the well- 



know n degeneracy between metallicity and age ( Wort hey 



1994) makes such conclusions very uncertain. 

For that reason, we have started a campaign to 
spectroscopically study the globular cluster system in 
NGC 3115. Our precision measurements of absorption 
line-strengths can be used to derive age and metallicity es- 
timates directly from the comparison with new stellar pop- 
ulation models. Unlike photometric methods, with spec- 
troscopy we are also able to explore element abundance 
ratios for the GCs. We compare our results with other very 
recently obtained spectroscopic samples of GCs in early- 
type galaxies: in the giant Fornax elliptical NGC 1399 



(Forbes et al. 2001), the Sa/Sb galaxy M81 (Schroder 



et al. ^002), the SBO galaxy NGC 1023 (Larsen fc Brodie 



2002) and the Sombrero galaxy NGC 4594 ( [Larsen et al, 
2002) 



The paper is organized as follows. In Section |^, the 
observations and their reduction are discussed. Section ^ 
presents the colour distribution of our sample while in 
Section ^ the treatment of abundance ratios and new stel- 



lar population models are investigated. Our results on 
abundance ratios, age and metallicity distributions for 
GCs in NGC 3115, the Milky Way (hereafter MW) and 
M 31 are presented in Section ^ with a general discussion 
in Section |^. We present our conclusions in Section 0. 



2. Observations and data reduction 

The candidate GCs were selected from the HST/WFPC2 
investigation of Kundu fc Whitmore| ( |1998| ) who detected 
144 globular cluster candidates in the central region of 
NGC 3115. In order to keep integration times reasonably 
short only clusters with V <22 (the peak of the GC lumi- 
nosity function is at F = 22.37±0.05; Kundu & Whitmorc 
p^98| ) were selected while keeping a balance between red 
((F - /) < 1.06) and blue {{V - I) > 1.06) clusters. In 
order to utilize the full field-of-view of our multi-object 
spectrograph we supplemented this list with GC candi- 
da tes from a low re solution spectroscopy survey reported 
in 



Kavelaars ( 1998 ) and also placed some slits on promis- 
ing objects without prior information (see Table |l]). In 
total two masks with each 26 slitlets covering more than 
50 different GC candidates were manufactured. 

The observations were carried out 26/27 Feb 2001 with 
F0RS2 on VLT using the blue 600 1/mm grism and 1" 
wide slitlets giving a resolution of -4.8 A (FWHM) sam- 
pled at 1.2 Apixel"^. The seeing was generally FWHM< 
I'.'O. The total exposure time for the first mask was 12 440 s 
split up in six individual exposures of varying length. Due 
to bad weather the second mask was only exposed for 
2700 s which was not long enough to be of use for this 
study. The flux standard GD 71 was observed through a 
long-slit (2'.'5 width) to enable us to correct the continuum 
shape of the spectra. 

The standard data reduction procedures (bias subtrac- 
tion, flat-fielding, wavelength calibration, sky-subtraction 
and continuum correction to a relative flux scale) were 
performed with a combination of MIDAS and IRAF tasks. 
For each slit a 2-dimensional subsection of the CCD was 
extracted and then treated as a long-slit observation. The 
extraction region was defined by tracing the globular clus- 
ter spectrum along the wavelength direction and extract- 
ing the corresponding sections for the flat-field and arc- 
lamp observations. The cosmic rays were removed with 



the routine lacos_spec (van Dokkum 2001). The wave- 
length calibration was accurate to < 0.2 A. After sky- 
subtraction GCs were extracted from the individual ex- 



posures to achieve an optimal S/N (Hornc 1986). Finally, 
the spectra of the individual exposures were combined. 

We used the spectroscopic flux standard GD 71 to 
correct the continuum shape of our spectra to a rel- 
ative flux scale and applied a reddening correction of 
E{B-V) = 0.146 dSchlegel et aL||1998| ). In order to trans- 
form our observations onto the Lick/IDS system we con- 
volved our spectra with a wavelength-dependent Gaussian 
kernel (taking into account small variations of the instru- 
mental resolution with position on the chip) thereby re- 
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producing the Lick/IDS spectral resolution (Worthey & 
Ottavi 



mi 



1997) 



As a first analysis step we measured the recession ve- 
locity (with fxcor in IRAF) and S/N of the GC candi- 
dates. The first slit (Slitlet ID: 1) did not show any object 
which we ascribe to a field distortion towards the edges of 
the field-of-view. Basic information for the remaining 28 
globular cluster candidates is listed in Table |l]. The first 
column in this table is an object identifier, the second col- 
umn indicates from which source the object was selected 
from. The following columns list the J2000.0 coordinates 
and / photometry of |Kundu fc Whitmortj ( |1998[ ). The 



read-out noise of the CCD and the errors in the ve- 
locity determination. The Lick system is based on non- 
flux calibrated spectra, so one expects small offsets (e.g.. 



seventh column shows pseudo (V — I) colours derived from 
our spectroscopy (see Section^ for details). Columns eight 
and nine list the radial velocity and the mean S/N per 
pixel respectively. Finally, column ten indicates whether 
the GC candidate is a spectroscopically confirmed mem- 
ber of the NGC3115 system (see next paragraph). 

The radi al velocity of NGC 3 115 is listed as 7 20 ± 
5 kms-i by ^mith et al] ( |2000| ). [Kavelaarsl ( |1998| ) find 
a mean velocity of 620 kms~^ and a velocity dispersion 
of (7 = 177 kms~^ from low resolution spectroscopic ob- 
servations of 22 globular clusters. We consider all objects 
with radial velocities between 200 kms~^ and 1300 kms~^ 
as members of the NGC3115 system. Object 6 is a back- 
ground galaxy with a redshift of z = 0.39, while object 4 
is likely to be a galactic star. In total 26 out of 28 GC can- 
didates (93%) turned out to be objects which are dynam- 
ically associated with NGC3115. This low contamination 
fraction is to be expected for a sample of candidates pre- 
dominantly selected by morphological criteria from HST 
data. 

In order to allow a proper line-strength analysis we 
restrict our sample to GCs with a mean S/N > 12 per 
pixel0 yielding a final sample of 17 GCs in NGC3115. 
Below this signal-to-noise cut the line-strength errors are 
too large to allow a proper age-metallicity analysis. The 
mean S /N of our final sample is 27 per pixel (ranging from 
12 to 61). The GCs have a mean velocity of 600 kms~^ 
(Icr scatter: 215 kms~^). The spectra of these 17 GCs 
broadened to the Lick/IDS resolution are shown in Fig. |l|. 

We note that the globular cluster 14a (IIST27, 
S/N=ll) shows clear [OIII]A5007 emission which is how- 
ever redshifted by «190 kms^^ with respect to the reces- 
sion velocity of the globular cluster itself. Therefore the 
emission cannot be associated with the globular cluster 
itself, although the spatial extent of the emission is con- 
sistent with the size of the GC on the sky. 



Lick indices (for index definitions see Worthey & 



Ottavi mi 1997; Tragcr et al. 1998) were measured from 
the resolution and continuum corrected spectra covering 
a common observed wavelength range of 3670-5500 A. 
Uncertainties in the indices were derived by Monte-Carlo 
simulations which take into account the photon noise. 



Kuntschner 2000) in the zero-point of the line indices. 



^ The S/N is determined as an average over the observed 
wavelength range. The minimum S/N per resolution element 
is Ri25. 



Since no stars in common with the Lick stellar library 
were observed during our observing run we use the offsets 
established by Vazdeki^ (1999) for a large flux-calibrated 
sample of stars ('Jones library') to transform our index 
measurements onto the Lick system. Since typical metal- 
licities for GCs are between —2.0 < [Fe/H] < 0.0, we 
calculated the offset by averaging the values for [Fe/H] = 
-0.7 and -0.4 given by |Vazdekis| ( |1999| , see his Table 2). 
Note, that some indices (e.g., Mg2) show strong evidence 
for a metallicity dependent offset, which can introduce sys- 
tematic offsets as a function of line-strength. Here in this 
paper we use mainly indices where the former problem is 
only a second order effect. The measured in dex v al ues a nd 
their associated errors are listed in Tables A.l & A. 2 in 
the Appendix. 

The positions of the GCs with respect to the parent 
galaxy are shown in Fig. |^. Due to the optimisation of 
the target efficiency and the observational constraints of 
the F0RS2 instrument the GCs lie in a narrow band '--^25" 
eastwards of the galaxy centre parallel to the major axis of 
NGC 3115. Both the red and blue sub-samples were evenly 
distributed across the CCD. 

In Figure |^ we assess the kinematics of all globular 
clusters which have recession velocities consistent with be- 
ing members of the NGC 3115 system. There is a strong 
signature of GC rotation along the major axis of the 
galaxy. We confirm here the weak signal of rotation for 
the red cluster population detected by Kavelaars (1998). 
In our sample there is no clear difference between red and 
blue clusters, both show an equally strong signal of rota- 
tion. 

Since our sample of GCs is dominated by clusters close 
to the location of the major axis (see Figure ||) we note 
that the sample is probably biased to globular clusters 
associated with the disk- formation of NGC 3115. 



3. Colour distribution 

9 out of 17 GCs for which we have good quality spec- 
tra and are kinematically associated with NGC 3115 are 
within the HST/WFPC2 image. In order to assign colours 
to the remaining 8 GCs we derived an empirical relation- 
ship between a "pseudo-colour" measured from our flux- 
calibrated spectra and the HST {V — I) colours. Of course 
our limited wavelength range does not allow us to mea- 
sure a real {V — I) colour but one can get a rough idea 
on the overall spectral shape. On the spectra we measured 
the mean counts in two regions (3900 - 4600 A and 4800 - 
5500 A) and plotted the difference versus the HST (V-I) 
colour (see Fig. ^). This procedure yields a linear relation- 
ship with a mean scatter of 0.037 mag in derived {V — 1) 
colour. Note that here we have used all available spectra of 
GC candidates with measured HST colours (see Table 
specifically including the ones with an S/N < 12 per pixel. 
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Rest wavelength [A] 

Fig. 1. Rest-frame spectra of globular clusters corrected to the Lick/IDS resolution and sorted by increasing strength 
of our mean metallicity indicator [MgFe] (from top to bottom, see also Section ^) . Note the change in overall continuum 
shape (metal poorer globular clusters are bluer) and also the change in absorption strength of the Balmer lines (HJ 
at 4101 A, H7 at 4340 A and H/3 at 4861 A). The slitlet ID for each spectrum is indicated at the right hand side of 
the plot (see also Table |l|). 
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Slitlet ID 


Origin'' 


RA J2000.0 DEC 


(HST) 


{V - I) (HST) 


(V - I) (spec)" 


cz 


S/N 


member 


1 


- 








No object in slit. 










2 




10''05'"06'. 


5 


-7°45'"01" 






1.183 


408zb37 


9 


V 


3 




10''05'"06'. 


1 


-7°44'"30" 






1.142 


439zbl7 


23 


V 


4 


- 


10''05"'09'. 


9 


-7°45'"Q1" 






0.883 


0±11 


72 


- 


5 


WHT09 


10''Q5'"11'. 


4 


-7°45'"Q1" 






0.980 


379zb21 


32 


V 


6 


- 


10''Q5'"Q9'. 


2 


-7°43'"52" 






- 


7=0.39" 


- 


- 


7 


HST08 


10''Q5™12'. 


6 


-7°44'"28" 


ZU.^O^ HI u.uuu 


1 1 70 _|_ n nno 


1.148 


291±16 


28 


V 


8 


HST03 


10''05'"12'. 


6 


-7°44'"09" 


1 Q 7Q/1 + n nriR 

ly. 1 m u.uuu 


±.uo± zc u.uuy 


1.029 


393±13 


61 


V 


9 


HST02 


10''Q5™13'. 


1 


-7°43"'59" 


1 Q 749 + 00^1 
ly. 1 ziz u.uuu 


1 (-177 _|_ n nn7 
±.U 1 ( zc u.uu / 


1.068 


494zbl6 


38 


V 


10 


HST21 


10''05'"13'. 





-7°43'"41" 


91 079 -1- n m 9 


n Q7Q -1- n m Si 
u.y ( y It u.uio 


1.050 


456zb37 


8 


V 


11 


HST39 


10''Q5'"14'. 


7 


-7°43'"56" 


91 f^'^Q -1- n ni /I 


1 1 QQ -1- n ni Si 


1.218 


702zb31 


8 


V 


12 


HST13 


10''Q5™14'. 


4 


-7°43'"35" 


ZU.Uy:^ ZIZ u.uuo 


n Q7Q -1- n ni 1 

u.yiy zc u.uii 


0.982 


582zb35 


19 


V 


13 


HST18 


10''05™15'. 


2 


-7°43'"34" 


on Q7n + n 01 n 

zu.y lU ZIZ u.uiu 


u.yoy zc u.uio 


1.015 


546zb28 


12 


V 


14a 


HST27 


10''Q5™16'. 


6 


-7°43"'29" 


91 9fifi + n ni 1 

Zr.ZUU ZIZ U.U-L_L 


1 1 fif^ -1- n ni 

i.iUO ZIZ U.UIO 


1.117 


321zb25 


11 


V 


14b 


HST12 


10''05'"17'. 


2 


-7°43'"20" 


on cioo -I- n nny 
zu.u^z ZIZ u.uu ^ 


1 1 CQ _|_ n nnq 
±.±oy ZIZ u.uuy 


1.138 


696zbl5 


23 


V 


15 


HST17 


10''Q5™16'. 


8 


-7°43'"09" 


on Q/ifi -1- n ni n 
zu.y^u ZIZ u.uiu 


n q7Q -1- n ni ^ 
u.y ( zt u.uio 


1.021 


739±23 


16 


V 


16 


HST32 


10''05™18'. 





-7°43'"15" 


01 /ifi7 -1- n ni 9 

AL.^Xjl ZIZ U.UIZ 


n qi Q -1- n ni 
u.yio zt u.uio 


0.878 


599zb54 


12 


V 


17 


HST09 


10''Q5™17'. 


2 


-7°42"'49" 


9n '^Q'^ -1- n 007 
zu.oyo zt u.uu/ 


1 1/17-1-0 nin 

l.l'^t/ zt U.UIU 


1.167 


813zbl7 


23 


V 


18 


HST31 


10''05'"18'. 


4 


-7°42'"45" 


21.392 ±0.012 


1.190 zt 0.016 


1.236 


210zb23 


9 


V 


19 


HST57 


10''05™17'. 


3 


-7°42'"08" 


21.921 ±0.018 


0.891 zb 0.028 


0.983 


736zblll 


5 


V 


20 


HST46 




6 


-7°41"'57" 


21.802 ±0.015 


0.919 zt 0.023 


0.946 


661zb55 


7 


V 


21 




10''05'"20'. 


1 


-7°42'"16" 






0.971 


846±42 


12 


V 


22 




10''^Q5'"19'. 





_7°41'"47" 






0.957 


670zb41 


7 


V 


23 


WHT15 


10'"05'"20'. 


2 


-7°41'"45" 






1.140 


682zbl5 


34 


V 


24 


WHT16 


10''Q5'"21'. 


2 


-7°41"'38" 






0.992 


1008zbl8 


42 


V 


25a 




10''05™23'. 


4 


-7°41'"59" 






0.832 


202zb42 


19 


V 


25b 




10''05'"23'. 


7 


-7°41"'55" 






0.890 


721zb30 


46 


V 


26a 




10''05'"22'. 


7 


-7°41'"13" 






1.068 


803zb24 


10 


V 


26b 


WHT17 


10''05™22'. 


9 


-7°41'"07" 






0.934 


769zb25 


19 


V 



Background galaxy where [On] emission is visible. 

HSTxx denotes globular c l usters selected from the HST study of Kundu & Whitmore (1998), WHTxx are objects from the 

(1998) while "-" indicates objects selected from our own VLT-FORS2 image in order to fill the 



Kavelaar! 



investigation of 
mask. 

Pseudo {V — I) colour derived from the globular cluster spectra. The mean error is 0.037 mag. See Section ^ for details. 



The derived {V — I) colours are listed in Table ^in column 
seven where we have assigned a mean error of 0.037 mag. 

The colour distribution of our sample of 17 GCs is 
shown in Fig. p| toge ther with the full HST sample by 
Kundu fc Whitmorcl (^998|). The HST sample shows a 
clear bimodal structure with peaks at {V — I) — 0.96 
and {V — 1) = 1.17. Taking the mean colour of the HST 
sample {{V — 1) = 1.06) as the dividing line between the 
red and blue clusters, our spectroscopic sample features 
eleven blue clusters and six red ones while spanning the 
range [V - I) = 0.87-1.18. 

4. Treatment of abundance ratios 

Theory of chemical enrichment predicts that stellar popu- 
lations created in a short burst of star formation show ele- 
vated magnesium-to-iron abundance ratios while extended 
periods of star formation lead to roughly solar abundance 



ratios (e.g., Tinsley 


1979 


• 


Worthey et al. 


1992; 


Worthey 


1998; Thomas et al. 


1999 


)• 


This is generally explained by 



either a delay in the production of Fe-peak elements with 
respect to a-elements due to the different timescales of 
SNIa and SNII, or by a star-burst with an initial mass 
function skewed to massive stars. For example, many ob- 
servations of bright elliptical galaxies in clusters show in- 



star formation (e.g., 


Peletier 


1989 


; IWorthey et al. 


1992 


Davies ct al. 


1993 


Kuntschnc 


I 2000). 



GCs represent simple stellar populations (hereafter 
SSP, i.e., a unique age and chemical composition) since 
all the stars of one GC are thought to be created in a 
single star formation event. But it is not known a pri- 
ori whether the stars of a given GC are formed out of 
Fe-deficient gas clouds, which have been only enriched by 
SN II producing little Fe, or out of an already well mixed 
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Fig. 2. Positions of spectroscopically confirmed globular 
clusters with respect to the centre of NGC3115 (trian- 
gles represent blue clusters and circles red clusters, see 
Section ^) . Approaching and receding globular clusters are 
indicated by open symbols and filled symbols respectively, 
while the symbol size indicates the extent of the velocity 
difference with respect to the mean velocity of the whole 
sample. The ellipse encloses half of the integrated light and 
indicates the position and orientation of the main galaxy 
itself (data from Michard fc Marchal 1994 ). North is up 



and east to the left. 

interstellar medium harbouring the products from both 
SN la (main producer of Fc) and SN II. For GC systems 
with a bimodal colour distribution like in NGC3115, all 
scenarios that have been proposed to explain the origin 
of the red metal-rich GCs start from the principle that 
the red population is formed in a separate star formation 



event (e.g., Ashman fc Zepi 1992 ; Forbes ct al. 1997 ). In a 
naive star-formation scenario, where the red clusters form 
from the well-mixed interstellar medium they should show 
solar abundance ratios. 

Most of the currently available stellar population mod- 
els which can be used to investigate the abundance ra- 
tios of extragalactic objects are based on stellar libraries 
from our own Galaxy. This has the disadvantage that, 
particularly at sub-solar metallicities, galactic disk stars 
show super-solar abundance ratios for many a-elements 
([Edvardsson et alj |1993|; |McWilliam| |99^. Therefore, 



without correction, the model predictions will be biased 



towards super-solar abundance ratios (e.g.. 


Borges et al. 


1995; Kuntschnci 


200C ; Thomas et al. 2002a 


)• 


ihomas et al. 


(2002a) (sec also Thomas et al. 


2002b) 



provide new models which take the stellar library biases 
into account and can predict line-strengths for solar abun- 
dance as well as non-solar abundance ratios over a large 
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Fig. 3. Assessment of rotation for NGC3115 globular 
clusters. The triangles and circles represent blue and red 
clusters, respectively. Filled symbols indicate the high S/N 
sample of 17 globular clusters for which we present a line- 
strength analysis in this paper. Panel a): The vertical axis 
shows the radial velocity of the globular clusters relative to 
the mean velocity of the sample. The horizontal axis is the 
position angle between a globular cluster and the galaxy 
centre, where a PA of degrees represents north and posi- 
tive is east of north. The sinusoid is a simple least-squares 
fit to the full sample with one iteration to reject outliers. 
Panel b): radial velocity is plotted against the projected 
major axis radius. Positive numbers are towards north- 
east. The small dots represent the stellar rotation curve 
of NGC 3115 along the major axis (Capaccioli et al. 1993). 



metalhcity range (—2.25 < [Fe/H] < 0.35). Since our pa- 
per is among the first that make use of these new mod- 
els for studying GC spectra, we first explore systemati- 
cally how the three parameters age, metalhcity and abun- 
dance ratio ([Mg/Fe]) affect the absorption line-strengths 
of SSPs. For this purpose we plot in Fig. || the model 
predictions for the line-strengths of the often used metal- 
hcity indicators Mg 6, (Fe)^ fc [MgFe]|^ and the age sensi- 
tive Balmer line H/3 as a function of metalhcity, age and 
[Mg/Fe]. 

The efl[ects of non-solar abundance ratios (thick lines) 
with respect to solar ratio model predictions (thin lines) 
can be clearly seen for the metalhcity sensitive indices 
Mg& and (Fe). At a given age and all metallicities 
(Fig. ||a), the Mg& index is predicted to be stronger 
for super-solar [Mg/Fe]. The difference between solar ra- 
tio and non-solar ratio models increases with increasing 
metalhcity. T his behaviour is rev ersed for the (Fe) in- 
dex (see also Trager et al. 2000| ). A similar effect can 
be observed when we fix [Fe/H] but vary age (Fig. ||b). 
Remarkably, the [MgFe] index, the geometric mean of the 
Mg6 and (Fe) indices, does not show any significant de- 
pendence on abundance ratios, at least not within the cur- 
rent framework of the models. This has the unique advan- 
tage that [MgFe] can be used as an empirical mean metal- 



^ (Fe) = (Fe5270 + Fe5335)/2 , ( |Gonzalez| [l99^) 

^ [MgFe] = ^Mgb X (Fe5270 Fe5335)/2 , (|Gonzami993|) 
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Fig. 6. Dependence o f Mgb, [MgFe], (Fe) and H/3 on [Fe/H] and age for models with different [Mg/Fe] as predicted by 
Thomas et al. ( 2002a| ). The thick lines and open diamonds stand for [Mg/Fe] = +0.3 whereas the thin lines and plus 



symbols represent [Mg/Fe] = 0.0. The respective line-style and average measurement error for the different indices is 
given on the left hand side of panels (a) and (c). For the [MgFe] index the two model lines (for [Mg/Fe] = 0.0 and 0.3) 
overlap almost completely and therefore the index is almost independent of the Mg to Fe abundance ratios. 



licity indicator with negligible dependence on the abun- 
dance ratios (see also Kuntschner et al, 200l| ). Similarly, 
H/3 is hardly affected by [Mg/Fe] with a small increase in 
line-strength for larger [Mg/Fe] ratios (Fig. ^&d). 

However, the H/? index shows a noteworthy compli- 
cation at low metallicities ([Fe/H] < —1.35) and ages 
>8 Gyrs. Here the models predict a non-monotonic de- 
crease of the H/3 index with increasing age (see Fig. ||d. 



models for [Fe/H] = -2.25 & -1.35). Therefore, in this 
age and metallicity range a given measurement of H/3 and 
[MgFe] indices does not correspond to one unique age but 
can in fact be consistent with a range of ages. The ambigu- 
ity arises between a genuine measurement of the turn-off 
temperature by H/3 and the appearance of blue horizontal 
branch stars in old, metal poor stellar systems which will 
start to increase the H/3 index, mimicking a younger age 
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Fig. 4. The pseudo colour derived from the spectra versus 
{V — /)hst- The linear fit corresponds to {V — /)hst = 
-1.08(±0.04) X (spec colour) + 0.95(±0.01). For details 
see Section \ 



Fig. 5. Histo gram of {V — I) colours of the cluster can- 



(Lee et al. 


200C; 


Maraston & Thomas 2000; 


Beasley et al. 


2002b 


) . The overall result is a "crossing" of iso-age lines at 



ages larger than 8 Gyr and low metallicities. In Section p. 2| 
where we will determine our best age estimates, the above 
effect will be taken into account. 



5. Results 

We first present our results on the abundance ratios of Mg 
to Fe of GCs in NGC3115 and compare these with the 



Lick/IDS observations (Frager et al 



1998) of GCs in our 



own galaxy and M31. In Section 5.2 we present our best 
estimates of the ages and metallicities of NGC3115 GCs, 
while we compare photometric and spectroscopic metal- 



licity estimates in Section 5.3 



5.1. A l^ndancc ratios 



In Fig. |7|, we show a diagram of Mgb vs (Fe) index. 
We plot our sample of GCs separated by {V — I) colour 
(see Fig. ^) where red and blue clusters are shown as 
open circles and open triangles respectively. Overplotted 
are the model predictions by Thomas et al. (2002a) for 
[Mg/Fe] = 0.0, 0.3, 0.5, age = 3, 5, 8, 12 Gyr and [Fe/H] = 
-2.25, -1.35, -0.33, 0.0, +0.35. In this diagram the effects 
of age and metallicity are almost completely degener- 
ate and the sensitivity towards abundance ratios is max- 
imised. The models predict a significant change in line- 
strength for varying abundance ratios at large metallic- 
ities/ages (i.e. large Mg6 & (Fe) line-strength) while at 
low metallicities the dependence is smaller. Therefore, at 



didates from Kundu fc Whitmorc ( 1998 ). Overplotted as 
a filled histogram is our spectroscopic sample of 17 con- 
firmed globular clusters in NGC 3115. The vertical dashed 
line indicates the dividing line between red and blue clus- 
ters at (V -I) = 1.06. 

a given error in line-strength, abundance ratios are deter- 
mined more accurately at larger metallicities/ages. 

We find that for GCs in NGC 3115 the abundance ra- 
tios vary from roughly solar to about [Mg/Fe] ~ -1-0.3 
(with a maximum of [Mg/Fe] ~ -1-0.5, see Fig. ^). In our 
small sample there is no evidence for a trend with colour. 
In particular, there are GCs with solar and super-solar 
abundance ratios in both the blue and red GC popula- 
tions. There is even some weak evidence for a bimodal 
abundance ratio distribution with peaks at [Mg/Fe] « 0.0 
and « 0.3 (see Fig. ^). This tentative evidence needs to 
be confirmed with a larger sample of GCs. 

The measurements for the centre of NGC 3115 itself 

and 



1998) 



(large filled square, data from Trager et al 
the radial g radient (small filled squares, data from Fishei 
et al. |1996[ ) along the major axis up to 40" is compatible 



with a model of [Mg/Fe] « 0.0 at high metallicity/age. 
In Fig.s 0b and c we plot the Lick/IDS observations of 
GCs in the Milky Way and M31 respectively (data from 
Trager et ar||1998| ). While virtually aU MW GCs are con- 
sistent with [Mg/Fe] ~ -1-0.3, similar to large elliptical 
galaxies, we find a range in [Mg/Fe] for the GCs in M31. 
The overall distribution of the [Mg/Fe] ratios for GCs in 
M31 is similar to the one we find in NGC 3115. The av- 
erage value of [Mg/Fe] = +0.3 we find for the MW GCs 
compares well with high resolution studies of individual 



stars in MW GCs (e.g., |Lee fc Carney||2002D . 

GCs with super-solar [Mg/Fe] ratios were previously 
found in other nearby galaxies but without our quan- 
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Fig. 7. Probing the [Mg/Fe] ratios of globular clusters in a Mg6 vs (Fe) diagram. Panel a): Our sample of globular 
clusters in NGC3115 is shown as open triangle s (blue cluste r s) and open circles (red clusters). The large filled square 
represents the centre of NGC3115 taken from Trager et al. (1998) and the small filled squares represent the data of 
Fisher et al. (1996) which cover radii up to 40" along the major axis. Panel b): Milky Way globular clusters observed 
with the Lick/IDS instrumentation. Overplotted as stars are early-type galaxies in the Fornax cluster with central 
velocity dispersion > 75 kms^^ ( [Kuntschner 2000). Panel c): Globular clusters in M31 observed with the Lick/IDS 
instrumentation. The filled triangle represents the centre of M31 taken from [Trager et al. (199S). Overplotted in all 
panels are models by Thomas et al. ( ^002a ) with abundance ratios of [Mg/Fe]=0.0., 0.3, 0.5 as indicated in the left 
panel. The models span a range in age (3-12 Gyr) and metalhcity ([Fc/H] — —2.35 to + 0.3). 



in NGC 1399. Using Mg and TiO features, Larsen et al 
(2002) find a mean [a/Fe] of +0.4 for both metal-poor and 



metal-rich GCs in the Sombrero galaxy. 
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Fig. 8. Histogram of abundance ratios [Mg/Fe] for 
NGC 3115 globular clusters. The abundance ratios were 
determined from Fig. Ma. 



titative accuracy. For example, Forbes et al, (2001) at- 



tribute super-solar abundance ratios to 4 out of 10 GCs 



5.2. Age and metallicity 

In the previous section we were able to determine the 
abundance ratios of GCs without knowing the age and 
metallicity since the latter parameters are almost com- 
pletely degenerate in a Mg6 vs (Fe) diagnostic diagram. 
However, our earlier discussion of the model systematics 
shows that we need to take the abundance ratios into ac- 
count in order to estimate the age and metallicity of the 



GCs (see also e.g., Trager et al, 

pooi| ). 



2000 



Kuntschner et al 



Principle age sensitive lines within our observed wave- 
length range are the Balmer lines H/3, H7, and US. For H7 
and H(5, the dependence on a-element to Fe ratio is yet un- 
known. H/3 is only marginally sensitive to abundance ratio 
variations, at least in comparison to our average observa- 
tional error. To further minimise the influence of abun- 
dance ratios, we employ as metallicity indicator [MgFe], 
that also shows no significant [Mg/Fe] dependence (see 
Fig. Within the accuracy of our data sample, an H/3 
vs [MgFe] diagram can therefore be used to estimate the 
ages and metallicities of our NGC 3115 GCs without being 
significantly affected by abundance ratios. 
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In Fig. H we show diagrams of [MgFe] versus the three 
Balmer hnes for our sample of NGC 3115 GCs (left panels) 
and the respective data for GCs in the Milky Way and 
M31 from the Lick/IDS observations (middle and right 
panels; the index values for H/3, Mg6 and (Fe) were taken 
from Trager et al. ( 1998| ); the higher order Balmer lines 
of the Lick/IDS observations are presented in Table A. 3 



The Lick/IDS observations of MW GCs also show a 
significant number of objects below the model predictions. 
We note that there are no systematic observational off- 
sets to be expected since the data was taken with the 
original Lick/IDS system. We speculate that observations 
of these GCs may be contaminated by fore/back-ground 
stars. New spectroscopic observations of MW GCs (Puzia 



abundance ratio models by Thomas et al. ( 2002aD and 
Maraston (2002, in preparation) for metallicities [Fe/H] = 
-2.25, -1.35, -0.33, 0.00, 0.35 (dashed lines, left to right) 
and ages 3, 5, 8 and 12 Gyr (sohd lines, top to bottom). 

We first caution that a direct comparison between 
models and our data to derive absolute ages and metallic- 
ities can be dangerous due to possible systematic calibra- 
tion errors. However, we estimate that systematic observa- 
tional errors are smaller than 0.1 A for the indices shown 
here and emphasize that relative comparisons within one 
sample will be significant. 

The observed H/3 values for NGC 3115 GCs show 
a large spread with respect to the model predictions. 
However, the data points which are well below the model 
predictions are the ones with the largest errors. Most 
of the well determined data points are close to the re- 
gion of a 12 Gyr SSP model. Our [MgFe] measurements 
show that there is a clear distinction in metallicity be- 
tween blue (open triangles) and red clusters (open circles) , 
with the red clusters being more metal rich (a more de- 
tailed analysis of the metallicity distribution is presented 
in Section 5.3). 

Since only about half of our data points have 
small enough error bars to be useful for an individ- 
ual age/metallicity evaluation we also calculate the er- 
ror weighted mean of the blue and red clusters, respec- 
tively. These average values (filled symbols in Fig. ^, 
left panels) give for the metal poor (blue) population an 
age of 12.0(1:^:^) Gyr and [Fe/H] = -1.05(±0.09), while 
the metal rich (red) population has an estimated age of 
lO-Sdli) Gyr and [Fe/H] = -0.26(±0.05). The errors on 
the age and metallicity are quoted as la errors on the 
mean values. 

We note in Section ^ that for metallicities [Fe/H] < 
— 1.35 and an age larger than 8 Gyr the strength of the H/3 
and [MgFe] indices is not uniquely connected to one age 
anymore. In fact there is a "crossing" of iso-age curves. For 
clarity we do not plot iso-age lines for ages greater than 
12 Gyr in Fig. ||, but this effect has been taken into account 
when deriving the errors on our best age and metallicity 
estimates. 

From the H/3 vs [MgFe] diagram we conclude that 
within our observational errors the two populations of GCs 
in NGC 3115 have the same age of 11-12 Gyr (assum- 
ing the calibration of models and data is accurate). The 
observed indices for the integrated light in the centre of 



& A.4| in the Appendix). Ovcrplottcd in Fig. 9 are solar- et al. 2002a) support this hypothesis since the new obser 



vations do not show such low H/3 values. Consistent with 
recent age estimates from the resolved stellar populations 



of M W GCs (e.g., [Rosenberg et al.||l999| [Salaris fc Weis: 
poo^ ) we do not find evidence for clusters younger than 
«8 Gyr. 

The Lick/IDS observations for M31 GCs show a rel- 
atively small scatter close to a 12 Gyr model predic- 
tion, with only three, metal-rich clusters showing evidence 
of a younger age. We note however, that for metallici- 
ties [Fe/H] < —1.35 the models seem to systematically 
over-estimate the H/3 absorption strength (or alternatively 
over-estimate the [MgFe] absorption strength). 

In the next paragraphs we will present our measure- 
ments of the higher order Balmer lines Hjf and HSp- 
We emphasize here that while these indices can be mea- 
sured with a higher precision than H/3, it is currently 
unknown how these indices depend on abundance ratios. 
Furthermore the absolute calibration of these indices has 
not yet been investigated in as much detail as the H/3 in- 
dex. 

The distributions of H-fp and HSf vs [MgFe] are nar- 
rower compared to H/3 vs [MgFe] and mostly encom- 
passed by the model grid. The error weighted means for 
NGC 3115 GCs indicate an age of w7 and «5 Gyr for 
the blue and red clusters, respectively. These average ages 
are substantially lower than what we inferred from the 
H/3 vs [MgFe] diagram. We note that our observations of 
H7i? and HSp for GCs in NGC 3115 agree well with the 
Lick/IDS observations of M31 and therefore we conclude 
that the calibration of the models is not consistent be- 
tween H/3 and the higher order Balmer lines. Despite this 
absolute calibration problem we find a good agreement in 
a relative sense between H/3, H'Jf and HSp- Therefore, at 
least to first order, we can say that the higher Balmer in- 
dices are not significantly affected by abundance ratios in 
the metallicity range probed by our data. 

Comparing the distributions for MW and M31 GCs, 
we find that the MW one is broader and offset towards 
smaller Hjf and HSp absorption consistent with the H/3 
measurements. We ascribe this to contaminated observa- 
tions for MW GCs (see above). The Lick/IDS data (partic- 
ularly the H/3 vs [MgFe] diagram), suggest that perhaps 
~3 metal rich M31 GCs have younger ages (3-5 Gyr). 
Alternatively, one could account for the strong H/3 absorp- 
tion in these metal rich clusters if the H/3 index is signifi- 
cantly influenced by an extended blue horizontal branch in 



NGC 3115 (taken from Trager et al. 199^, and shown as 
filled square in Fig. |9|) are consistent with a luminosity 
weighted age of Ril2 Gyr. 



an otherwise old, metal rich stellar population. Maraston 



& Thomas (200C) show that this effect can play a role in 



metal poor stellar populations, however, to date there is 
scarce observational evidence for the existence of a popu- 
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Fig. 9. Age and metallicity diagnostic diagrams using as metallicity indicator [MgFe] and as age indicator H/3, Hjf, 
and HSf- H/3 and [MgFe] are not significantly influenced by abundance ratios, for H'Jf, and HSf the behaviour is 
unknown. Our sample of globular clusters in NGC 3115 is shown in the left panels as open triangles (blue clusters) and 
open circles (red clusters). The error weighted means of the blue and red cluste rs are shown as filled symbols. The filled 
square in the left panels represents the centre of NGC 3115 Trager et al. ( 199g ). The middle and right panels show the 
Lick/IDS data for Milky Way MW and M31 globular clusters, respectively. Note, that for the MW globular cluster 
NGC 2158 there are n o Hjf and HSf index measurements available. Overplotted are solar-abundance SSP models by 
frhomas et al] ( ^002a| ) and Maraston (2002, in preparation) for metallicities [Fe/H] = -2.25,-1.35,-0.33,0.00,0.35 
(dashed lines, left to right) and ages 3, 5, 8, and 12 Gyr (solid lines, top to bottom). 
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lous extended blue horizontal branch in metal rich clusters. 
Rich et al. (1997) detected a blue horizontal branch in two 



metal rich MW GCs and Ferraro et al.| ( ^001 ) detected 



UV-excess stars in the core of 47 Tuc (see also Moehler 
et al.| ^OOOD 



Few spectroscopic observations of GCs in early-type 
galaxies with sufhcient S/N to investigate these effects 
have been published. Forbes et al. (2001) find that most 
of their 10 GCs in the giant elliptical NGC 1399 are old 
and compatible with a model age of 11 Gyr (using models 
by Maraston 2002, in preparation). Only two GCs dis- 
play such large H/3 values that these have either a very 
young age of ~ 2 Gyr or are "contaminated" by a signifi- 
cant blue horizontal branch population which causes large 
H/3 absorption. The aut hors prefer the first interpretation. 
Larsen fc Brodie (2002) present spectra of 14 GCs in the 
Sombrero Galaxy (NGC 4594). Their analysis of the co- 
added spectra of metal-poor and metal-rich GCs leads to 
age estimates between 10-15 Gyr. The m ajority (11 out 



of 14 GCs) of the spectroscopic sample of Schroder et al 
( 2002| ) of M81 GCs is c ompatible with old ages (using 



models by |Worthey| [19941) . There is only one outlier with 
a very high H/? line strength. 

In summary we conclude from our best calibrated dia- 
gram of H/3 vs [MgFe] that the majority of our sample of 
GCs in NGC 3115, regardless of their metallicity, are con- 
sistent with an age of «12 Gyr. Only one, metal rich clus- 
ter (Slitlet ID: 7) shows a combination of II/3 and [MgFe] 
absorption strength which indicates an age lower than 
8 Gyr. The higher order Balmer lines indicate a narrow 
distribution in age, with a hint of the metal rich clusters 
being younger by «2 Gyr. The unknown dependence of 
the higher order Balmer lines on abundance ratios makes 
this age difference highly speculative. The absolute ages 
indicated by the higher order Balmer lines are lower com- 
pared to the H/3 index. We ascribe this age difference to 
an inaccurate calibration of the higher order Balmer lines 
in the current stellar population models. The Lick/IDS 
samples of MW and M 31 GCs also show old stellar pop- 
ulations; only ~3 GCs in M31 show tentative evidence of 
younger stars. 

5.3. Photometric versus spectroscopic metallicity 
estimates 

In this section we compare our spectroscopic metallicity 
estimates with photometric methods and also investigate 
the general distribution of metallicities. For this purpose 
we assume an average age of the GCs in NGC 3115 of 
12 Gyr which is consistent with our findings in the previ- 
ous section. 



Fig. |10| a shows the purely empirical relation between 
{V — I) colour and our mean metallicity indicator [MgFe] . 
There is a tight relation over the observed parameter 
space. Overplotted as solid line are model predictions 



by Maraston (2002, in preparation) and Thomas ct al 



(2002a) for a constant age of 12 Gyr, which is in excellent 



agreement with our data. We note that the model pre- 
dictions for colours do not include the effects of non-solar 
abundance ratios. 

In order to convert the colours into metallicity esti- 
mates several authors have derived linear conversion for- 
mulae based on observations for MW GCs. For exam- 
ple, [Kundu fc Whitmorc] ( |1998| ) conclude that [Fe/H] 
= — 5.89 + 4.72(T^— /) is a good linear approximation. One 
can also use the predictions of stellar population models 
(Maraston 2002, in preparation) to predict the relation 
between {V — I) colour and metalhcity [Fe/H]. A com- 
parison of the empirical and synthetic calibration (12 Gyr 
model) is shown in Fig. ^0|b. 

Overall, the agreement is acceptable, although there 
are significant differences. Specifically at the low metal- 
icity end the models predict a shallower trend than the 
empirical relation. In order to stay consistent with the 
Kundu fc Whitmorc ( 1998 ) paper we use their relation to 
convert {V — I) colour to metallicity (see also Table Q). 
Furthermore we determine metallicity estimates from our 
spectra by using the [ MgFe] index in conjimction with the 
model predictions by Thomas et al. (2002a) and assuming 
a constant age of 12 Gyr. The comparison between pho- 
tometric and spectroscopic metallicity estimates is shown 
in Fig. |lO| c. 

Wc find a good linear relation between both meth- 
ods. The best fitting linear relation including the observa- 
tional errors is [Fe/}l]^^^^^ = -0.26(±0.05)-t-0.95(±0.08) x 
[Fe/II]gpj,^ with a probability of 0.30. The systematic 
offset of approximately [Fe/H] = —0.26 in the sense that 
the spectroscopic metallicity measurements are larger is 
consistent with the difference between model predictions 
and empirical calibration of the conversion formulae be- 
tween colour and [Fe/H] as shown in Fig. |lO|b. The pre- 
dicted non-linearity of {V — /) colour as function of metal- 
licity below [Fe/H] = — 1.5 cannot be tested since our data 
do not really cover this range. 

In summary we confirm with our accurate spectro- 
scopic observations that the bimodal colour distribution 
seen in NGC 3115 GCs is dominated by a metallicity ef- 
fect rather than by an age difference. Furthermore, both 
{V — I) colour peaks do show a substantial spread in 
metallicity. We conclude that in the metallicity range 
— 1.5 < [Fe/H] < 0.0 and in absence of young GCs, the 
{V — I) colour is indeed a good indicator for metallicity. 
We note that for metallicities below [Fe/H] = —1.5 this 
may not be the case anymore. 

6. Discussion 

Before we discuss the implications of our results it needs 
to be noted that our sample of GCs in NGC 3115 is dom- 
inated by the bright end of the luminosity function. Our 
faintest object (Slitlet ID: 17) which we could use for the 
line-strength analysis is 0.9 mag brighter than the peak 
of the GC luminosity function. Furthermore, since this 
project was only a pilot study we observed only a small 
number of clusters over a limited region across NGC 3115. 
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Fig. 10. Comparison of photometric and spectroscopic metallicity estimates. In the left panel the [MgFe] index is 
plotted against {V — I) colours (including the pseudo {V — I) colours derived from the spectra). Overpl otted as 
solid line is t he prediction for a 12 Gyr stellar population model by Maraston (2002, in preparation) and Thomas 
et alj_^ 2002a). In the middle panel we show a comparison between an empirical calibration of {V — I) colour against 



metallicity (dashed line; Kimdu & Whitmore 1998) and predictions from stellar population models (solid line, 12 Gyr; 
Maraston 2002, in preparation). In the right panel, photometric metallicities are calculated according to Kundu & 
Whitm|or(^ (|l998[), wh ile spectroscopic metallicities are derived from the [MgFe] index in comparison to models of 
Thomas et al. (2002a) assuming a constant age of 12 Gyr. The data point at [Fe/Hjgpoc [MgFe] = —2.25 represents 
a globular cluster with line-strength weaker than the currently available model predictions and is therefore only an 
upper limit. The lines at the right hand side of the plot show the peak positions of the bimodal distributions of GCs 
observed by Kundu & Whitmore (1998, solid lines) and RGB stars observed by Elson (1997, dotted lines). 



Although our data can give clear insights into the globular 
cluster system of NGC3115 it is not necessarily a repre- 
sentative sample. The comparison with the Lick/IDS ob- 
servations of MW and M31 GCs needs to be viewed with 
caution since these samples cover a larger range in abso- 
lute GC magnitude. 



Using their HST V, I photometry Kundu & Whitmore 
( |1998| ) develop the following formation scenario for 
NGG 3115. The blue, metal poor clusters are formed with 
the halo/bulge component of the galaxy very early on. 
Then about 4 ± 3 Gyr later an unequal mass, gas-rich 
merger event forms the disk component and the associated 
red, metal rich clusters. The authors point out that the 
spatial distribution of the blue and red clusters are con- 
sistent with the halo/bulge and disk components respec- 
tively. Furthermore, there is evidence from optical imaging 
( ^ilva et al. 



, that the disk is bluer and therefore per- 
haps younger than the halo/bulge component, consistent 
with the above outlined scenario. 

In our spectroscopic study of the GCs in NGG 3115 
we find that the clusters are consistent with being co- 
eval at about 11-12 Gyr. There is perhaps a weak hint 
of the red clusters being younger, but by no more than 
~2 Gyr. This on its own would not rule out the scenario 



by Kundu & Whitmore (1998), however our estimates of 



the abundance ratios show that at least the population 
of red, metal rich clusters is not homogeneous. For these 
objects we find a range in abundance ratios from solar to 
about [Mg/Fe] ~ 0.5. In our small sample of limited spa- 
tial extend we do not find any clear trends of the chemical 



composition of GGs with the kinematics (see Figure g). 
The relative velocities of most of the clusters in our sam- 
ple are consistent with rotation. A larger sample of more 
complete spatial coverage is needed to establish possible 
trends between abundance ratios and e.g., age, metallicity, 
spatial position and kinematics. 

What we can however say is that the metal-rich clus- 
ters with solar abundance ratios must have been made 
out of well mixed material which incorporates the prod- 
ucts of both SN II (the main producer of alpha elements) 
and SNIa (the main producer of Fe-peak elements). Since 
SNIa are somewhat delayed compared to SNII the solar 
abundance ratio clusters must have formed after the initial 
star burst in NGC 3115. There are many possible scenar- 
ios to explain the observed abundance ratio distribution, 
but it is hard to fit them into a simple picture of only two 
distinct formation events which create the red and blue 



globul ar cluster sub-populations (see also [Beasley et al 
^002a| ). 



One scenario which we would like to put forward for 
further discussion is the following. The metal poor (blue), 
non-solar abundance ratio clusters are associated with the 
halo formation, and the metal-rich (red), non-solar abun- 
dance ratio clusters are formed together with the bulge as 



was similarly proposed for the Milky Way by Carney et al 
( |1990D and |Wyse fc Gilmorc| ( |1992D . The metal rich, solar 



ratio clusters are then formed with the disk of this lentic- 
ular galaxy '^1-2 Gyr after the initial star-burst perhaps 
in connection with a merger. This scenario would then re- 
quire the disk to have also close to solar abundance ratios, 
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which can be tested by observations of the integrated light. 
Furthermore, if this scenario is correct the spatial distribu- 
tion and kinematics of disk GCs will be distinct from the 
population of halo and bulge GCs in NGC3115. Future, 
larger samples of NGC 3115 GCs will be very valuable to 
explore the connection between disk formation and metal 
rich globular cluster formation. 

More spectroscopic observations of GCs in nearby 
galaxies of various types will be very valuable to im- 
prove our stellar population models and learn more about 
the early star-formation epochs in early-type galaxies. 
However, the mismatch of the models and some observed 
indices demonstrates that we are also in urgent need for 
a new, high-quality flux-calibrated spectral library in or- 
der to exclude simple observational offsets. Only then can 
we make good progress with improvements on the input 
physics of stellar population models and their application 
to extragalactic objects. 

7. Conclusions 

We present new, accurate measurements of absorption 
line-strength indices of 17 globular clusters (hereafter 
GCs) in the nearby SO galaxy NGC 3115. Our objects span 
a range in colour so that the biniodal (V — /) colour dis- 
tribution is well sampled. 



A critical comparison with Lick/IDS data ( Trager 

let al.||l998| ) of GCs in M31 and the Milky Way (hereafter 
MW) is presented. The Lick/IDS measurements of the 
H7A,F and H(5a,_f indices for MW and M 31 GCs are pre- 
sented for the first time in this paper. The data are anal- 
ysed with new stellar population models (Thomas et al 
2002a) which are able to predict line-strength not only as 
a function of age and metallicity but also as a function 
of abundance ratio. Specifically, abundance ratio biases in 
the stellar library, which is an essential ingredient to the 
model predictions, have been taken into account for the 
first time. 

Our main results are listed in the following: 

— The GCs in NGC 3115 show a range of abundance ra- 
tios as estimated by the strength of Mg and Fe lines. 
Specifically we find for both red and blue clusters solar 
as well as super-solar values (up to [Mg/Fe] ~ -1-0.5). 
Lick/IDS data of M 31 GCs show a similar distribution 
while MW GCs are consistent with a constant value of 
[Mg/Fe] ~ 0.3. The latter is in agreement with recent 
studies of the resolved stellar populations in MW GCs. 

— Our analysis of II/3, Il7i? and HSp vs [MgFe] 
age/metallicity diagnostic diagrams shows that the red 
and blue GC populations are coeval within ?»2 Gyr. 
Our best calibrated diagram (II/3, [MgFe]) indicates 
a mean age of 11 — 12 Gyr. However, the higher or- 
der Balmer lines, although confirming the similar ages 
of blue and red clusters, indicate younger absolute 
ages («6 Gyr). We ascribe these younger age esti- 
mates to inaccurate model calibrations of the higher 
order Balmer lines. Evidence for the existence of young 



(< 5 Gyr) GCs in the studied galaxies is scarce. Only 
one cluster in NGC 3115 and perhaps ~3 clusters in 
M31 show a combination of Balmer and metal ab- 
sorption strength which is consistent with such young 
ages. The strong H/3 absorption can be alternatively 
explained if these clusters have old stellar populations 
with a populous extended horizontal branch. 

— We present a comparison of photometric and spec- 
troscopic metallicity determinations and find a good 
linear relation in the metallicity range probed by our 
sample of NGC3115 clusters (-1.5 < [Fe/H] < 0.0). 
The photometric estimates are systematically lower 
(~ —0.26) in comparison with our spectroscopic mea- 
surements. We note that our observations clearly show 
that each colour peak has a significant spread in metal- 
licity rather than being consistent with a narrow dis- 
tribution. 

— The existence of solar as well as elevated Mg-to-Fe ra- 
tios at a given metallicity for GCs in NGC 3115 in- 
dicates that a simple scenario of two distinct star- 
formation episodes is not sufficient to explain the for- 
mation of this galaxy. Probably a realistic model needs 
to incorporate more than two distinct star-formation 
events. 

— We detect a clear signal of rotation in our sample of 
GCs independent of their metallicities. 

Larger samples of high signal-to-noise spectra of GCs 
in nearby galaxies are required to increase our knowledge 
of how these galaxies formed. Particularly the measure- 
ment of abundance ratios in extragalactic GCs is now 
possible and will deliver important new constraints for 
galaxy formation models which are not accessible with 
broad-band colours alone. 
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Appendix A: Tables of line-strength 
measurements 

A.l. Line-strength measurements for globular clusters 
in NGC3115 



Table AJ & A^ present the line-strength measurements 
for globular clusters in NGC3115. The spectra have been 
broadened to the Lick/IDS resolution prior to index mea- 
surements. Small offsets have been made to match the Lick 
system as closely as possible (see Section |]for detai ls). For 
a summary of the definition of Lick/IDS indices see Trager 
let al.| (|1998| ). 



A. 2. Lick/IDS observations of H'yA,F ^nd H6a,f 
indices for Milky Way and M 31 globular clusters 

Since the higher order Balmer indices of Milky Way and 
M 31 globular clusters as measured in the Lick/IDS system 
have not yet been published, we list the measurements of 
H7A,_F and RSa.f in T ables |0| fc |a!^ . The Lick obser- 
vations are described in Burstein et al. (1984), names are 



explained there (see also Trager et al. 1998). M31 globular 
clusters were observed with the standard aperture of the 
Lick/IDS system (1'.'4 x 4'.'0). The MW globular clusters 
were observed with a long-slit of standard width that was 
raster-scanned on the sky to create a square aperture of 
size 66" X 66" (only for NGC 6624 the aperture is 45" x 
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Table A.l. Line-strength measurements for NGC3115 globular clusters 



Slitlet ID 




RSf 


CNi 


CN2 


Ca4227 


G4300 


RjA 


H7F 


Fe4383 


Ca4455 




[A] 

L J 


[A] 

L J 


[mag] 


[mag] 


[A] 

L J 


[A] 

L J 


[A] 

L J 


[A] 

L J 


[A] 

L J 


[A] 

L J 


3 


-1.41 


0.32 


0.049 


0.076 


0.95 


4.64 


-4.76 


-0.92 


3.74 


1.05 




0.54 


0.37 


0.015 


0.018 


0.26 


0.43 


0.51 


0.31 


0.59 


0.29 


5 


1.44 


1.81 


-0.014 


0.022 


0.54 


2.84 


-0.92 


1.04 


2.03 


0.40 




0.34 


0.23 


0.010 


0.012 


0.18 


0.30 


0.32 


0.18 


0.44 


0.21 


7 


-1.00 


0.29 


0.078 


0.120 


1.41 


4.45 


-4.60 


-0.53 


3.81 


0.67 




0.46 


0.32 


0.013 


0.015 


0.21 


0.35 


0.38 


0.23 


0.49 


0.22 


8 


0.49 


1.38 


0.006 


0.038 


0.83 


3.29 


-1.92 


0.54 


2.36 


0.56 




0.18 


0.12 


0.005 


0.006 


0.09 


0.15 


0.16 


0.10 


0.22 


0.11 


9 


0.68 


1.27 


0.015 


0.054 


0.79 


3.91 


-3.15 


-0.17 


3.09 


0.73 




0.28 


0.20 


0.008 


0.009 


0.14 


0.23 


0.27 


0.17 


0.35 


0.17 


12 


2.42 


2.87 


0.010 


0.067 


0.60 


2.38 


-0.39 


1.24 


2.53 


0.77 




0.52 


0.34 


0.016 


0.019 


0.27 


0.49 


0.51 


0.30 


0.71 


0.35 


13 


2.02 


1.95 


-0.011 


0.043 


0.94 


2.64 


-1.78 


1.20 


2.62 


0.04 




0.81 


0.57 


0.023 


0.028 


0.42 


0.72 


0.81 


0.49 


1.08 


0.53 


14b 


0.75 


1.23 


0.012 


0.035 


1.01 


3.93 


-4.56 


-0.58 


3.89 


1.26 




0.50 


0.36 


0.015 


0.018 


0.25 


0.40 


0.46 


0.28 


0.60 


0.28 


15 


2.54 


1.96 


-0.015 


0.020 


1.16 


3.29 


-1.72 


0.81 


3.00 


0.12 




0.64 


0.45 


0.019 


0.023 


0.34 


0.56 


0.64 


0.38 


0.82 


0.42 


16 


3.50 


2.95 


-0.063 


-0.058 


-0.21 


0.55 


0.96 


2.24 


2.07 


0.60 




0.84 


0.55 


0.025 


0.030 


0.45 


0.83 


0.78 


0.50 


1.11 


0.54 


17 


-0.27 


0.14 


0.056 


0.080 


1.07 


5.02 


-4.59 


-1.11 


3.53 


0.47 




0.54 


0.37 


0.014 


0.018 


0.24 


0.40 


0.48 


0.29 


0.57 


0.28 


21 


2.01 


2.63 


-0.048 


-0.034 


0.70 


3.34 


-0.17 


2.03 


1.85 


0.17 




0.85 


0.59 


0.024 


0.029 


0.43 


0.79 


0.79 


0.48 


1.11 


0.54 


23 


-1.79 


0.40 


0.069 


0.099 


0.90 


4.38 


-3.97 


-0.21 


4.77 


1.13 




0.37 


0.24 


0.010 


0.012 


0.16 


0.28 


0.31 


0.19 


0.37 


0.20 


24 


1.99 


2.06 


-0.008 


0.012 


0.83 


2.87 


-1.23 


0.93 


2.45 


0.32 




0.24 


0.17 


0.008 


0.009 


0.12 


0.22 


0.24 


0.14 


0.32 


0.15 


25a 


1.41 


1.53 


-0.049 


-0.041 


0.32 


2.66 


-0.31 


0.77 


0.46 


0.38 




0.54 


0.36 


0.015 


0.019 


0.28 


0.49 


0.51 


0.32 


0.71 


0.33 


25b 


2.33 


2.48 


-0.059 


-0.034 


0.49 


2.33 


0.30 


1.79 


1.27 


0.31 




0.22 


0.15 


0.006 


0.008 


0.11 


0.19 


0.21 


0.13 


0.31 


0.15 


26b 


2.46 


2.07 


-0.064 


-0.047 


0.54 


1.99 


1.35 


2.19 


1.33 


0.03 




0.53 


0.35 


0.016 


0.019 


0.28 


0.49 


0.49 


0.30 


0.71 


0.34 



Notes: Listed are line-strength measurements in the wavelength range 4000-4500 A for globular 
clusters in NGC3115 (see also Table |l|). The first line for each object shows the measured index 
values whereas the second line lists the associated errors. The unit for each index is given at the 
top in square brackets. 



60"). Raster scans have the same spectral resolution as 
standard slit width scans. For the definition of the higher 
order Balmer indices see Worthey & Ottaviani (1997). 
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Table A. 2. Lick/IDS line-strength measurements for NGC3115 globular clusters 



Slitlet ID 


Fe4531 


C24668 


H/3 


Fe5015 


Mgi 


Mg2 


Mgfo 


Fe5270 


Fe5335 


Fe5406 




L J 


L J 


L J 


L J 


fmael 

L &J 


[mae:l 


fAl 

L J 


L J 


L J 


L J 


3 


2.94 


3.88 


1.59 


5.24 


0.068 


0.216 


3.59 


2.30 


1.90 


1.46 




0.44 


0.65 


0.25 


0.56 


0.006 


0.007 


0.28 


0.28 


0.31 


0.23 


5 


1.87 


1.45 


2.12 


3.51 


0.038 


0.105 


1.86 


1.49 


1.05 


0.68 




0.34 


0.49 


0.19 


0.41 


0.005 


0.005 


0.20 


0.22 


0.25 


0.19 


7 


3.20 


4.07 


2.12 


5.90 


0.064 


0.196 


3.84 


2.03 


2.57 


1.74 




0.35 


0.53 


0.20 


0.43 


0.005 


0.006 


0.22 


0.24 


0.25 


0.20 


8 


2.20 


1.87 


2.18 


3.20 


0.025 


0.110 


2.04 


1.83 


1.73 


0.89 




0.17 


0.25 


0.10 


0.23 


0.002 


0.003 


0.10 


0.12 


0.13 


0.11 


9 


2.58 


2.89 


1.71 


4.13 


0.050 


0.160 


3.32 


1.97 


1.51 


0.83 




0.26 


0.39 


0.16 


0.34 


0.004 


0.004 


0.16 


0.18 


0.20 


0.16 


12 


0.70 


—0.29 


2.33 


3.53 


0.025 


0.118 


2.42 


1.88 


1.21 


1.21 




0.61 


0.86 


0.32 


0.69 


0.008 


0.010 


0.35 


0.37 


0.43 


0.34 


13 


1.99 


1.99 


1.38 


2.59 


0.014 


0.107 


2.12 


2.73 


0.90 


0.46 




0.87 


1.30 


0.49 


1.05 


0.011 


0.013 


0.52 


0.54 


0.63 


0.52 


14b 


3.26 


1.87 


1.38 


3.75 


0.052 


0.174 


2.90 


2.45 


2.69 


1.22 




0.43 


0.67 


0.26 


0.57 


0.006 


0.007 


0.27 


0.28 


0.31 


0.26 


15 


1.68 


—0.36 


1.31 


3.85 


0.035 


0.125 


2.25 


1.36 


1.87 


0.47 




0.66 


1.00 


0.37 


0.80 


0.009 


0.010 


0.39 


0.41 


0.48 


0.41 


16 


1.64 


-1.13 


2.65 


1.62 


0.015 


0.028 


0.56 


0.36 


0.30" 


-0.28 




0.92 


1.48 


0.56 


1.18 


0.013 


0.015 


0.57 


0.65 


1.12 


0.57 


17 


3.18 


3.47 


1.72 


4.80 


0.055 


0.199 


3.79 


2.24 


2.59 


1.19 




0.43 


0.66 


0.26 


0.52 


0.006 


0.007 


0.26 


0.28 


0.30 


0.25 


21 


0.66 


0.26 


1.39 


1.48 


0.009 


0.070 


1.07 


0.99 


1.13 


1.11 




0.92 


1.37 


0.53 


1.11 


0.012 


0.013 


0.53 


0.58 


0.63 


0.53 


23 


3.21 


2.36 


1.78 


4.31 


0.036 


0.166 


2.82 


2.68 


2.19 


1.65 




0.30 


0.45 


0.18 


0.36 


0.004 


0.005 


0.18 


0.18 


0.21 


0.16 


24 


2.01 


0.46 


1.96 


3.04 


0.023 


0.103 


1.71 


1.51 


1.31 


1.06 




0.25 


0.39 


0.15 


0.32 


0.003 


0.004 


0.16 


0.17 


0.19 


0.15 


25a 


0.12 


-0.86 


2.57 


0.97 


-0.004 


0.037 


1.44 


1.29 


0.78 


-0.33 




0.58 


0.87 


0.32 


0.81 


0.008 


0.009 


0.33 


0.38 
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0.66 
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0.32 



Notes: Listed are line-strength measurements in the wavelength range 4500-5450 A for globular 
clusters in NGC3115 (see also Table |l]). The first line for each object shows the measured index 
values whereas the second line lists the associated errors. The unit for each index is given at the 
top in square brackets. 

° For this globular cluster we measured a Fe5335 index of — 1.24(±0.75) A. This negative number 
is caused by a noisy section in the spectrum. By comparing this index with other Fe indices we 
established that a realistic value of the Fe5335 index for this globular cluster should be -1-0.3. To 
reflect the uncertainty in this procedure we increased the error by a factor of 1.5. 
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Table A. 3. Lick/IDS line-strength measurements for 
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Table A. 4. Line-strength measurements for globular clus- 
ters in M31 



Notes: We list here the Lick/IDS measure- 
ments of the indices H5a, H(5f, H7A and H7F 
for MW globular clusters. The Lick/IDS sys- 



t em i s described in detail in Trager et al 
( 1998 ). For the definit ion of the higher or- 



der Balmer indices see 



Worthey fc Ottaviani 



(1997). The first line for each object shows the 
measured index values whereas the second line 
lists the associated errors. The unit for each 
index is given at the top in square brackets. 
Three dots indicate where no index measure- 
ment is available. 
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Notes: We list here the Lick/IDS measurements 
of the indices RSa, RSf, H7A and H7F for glob- 
ular clusters in M31. The Lick/ID S syst em is 



described in detail in Trager et al. ( 1998 ). For 



the definition of the highe r order Balmer indices 



see Worthey fc Ottaviani (1997). The first line 



for each object shows the measured index values 
whereas the second line lists the associated er- 
rors. The unit for each index is given at the top 
in square brackets. Three dots indicate where no 
index measurement is available. 



